Design of coupling for targeting synchronization in time-delayed systems 
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We report a general method of designing delay coupling for targeting a desired synchronization 
state in delay dynamical systems. We are able to target synchronization, antisynchronization, lag-, 
anti-lag synchronization, amplitude death and generalized synchronization in mismatched oscilla- 
tors. We apply the theory for targeting a type of mixed synchronization where synchronization, 
antisynchronization and amplitude death coexist in different pairs of state variables of the coupled 
system. By this method, a scaling of the size of an attractor is possible in different synchronization 
regimes. We derive the stability conditions of synchronization based on Kravoskii-Lyapunov function 
theory and Hurwitz matrix criterion. We present numerical demonstration using the Mackey-Glass 
system and a delay Rossler system. 

PACS numbers: 05.45.Xt, 42.65.Sf, 05.45.Pq 



I. INTRODUCTION 

In conventional studies [1] of chaos synchronization, 
the coupling configuration is always known a priori, as 
examples, diffusive coupling, synaptic coupling. The sta- 
bility of synchronization is established above a critical 
coupling based on the negative measure of transverse 
Lyapunov exponent [2] . As a general condition, a master 
stability function [3] is used for determining the critical 
coupling strength for synchronization in identical oscilla- 
tors. Alternatively, an engineering approach to synchro- 
nization is recently attempted [4-9] for known systems 
where the coupling is unknown a priori. The coupling is 
rather designed that can ensure the asymptotic stability 
of the error dynamics of the coupled oscillators. So far a 
few specific issues were investigated such as engineering 
coherent patterns [5] , a desired amplitude dynamics with 
a target frequency [6] , synchronization in cellular systems 
[7] but each uses a separate approach to define the cou- 
pling. These studies were mainly restricted to non-delay 
systems and made in absence of coupling delay. How- 
ever, a design of coupling in the form of delay feedback 
[8] was proposed for effective control of synchronization 
and a desired desynchronization in non-delay system. A 
design of coupling was also proposed [9] for enhancing 
anticipation synchronization in presence of coupling de- 
lay in non-delay systems. A coupling delay [10] appears 
during the transmission of signals or information in many 
natural systems such as long axons of neurons [11] and, 
also in practical systems like communication channel [12]. 
This induces additional instability in the process of an 
emergent synchronization state. On the other hand, in 
physical systems like optoelectronic devices [13], even in 
electronic systems, time delay is an inherent property of 
the system. As example, consider a time-delayed system 
i.e., x = f(x(t),x(t — t)) where r is a delay time. This 
intrinsic time delay makes a system infinite dimensional 
and more complex. Synchronization in delay systems was 
thus given a special attention [10, 12,14-16] using delay 



diffusive coupling or a special type of delay coupling [16]. 
Stability of synchronization in delay systems was usually 
established using the Krasovskii-Lyapunov function the- 
ory [17] that showed dependence of system parameters 
and coupling strength. In all these studies, the coupling 
configuration was always pre-defined. An adaptive cou- 
pling was proposed recently for synchronization [18] in 
time-delayed optoelectronic delay feedback loops, where 
the coupling configuration was assumed known but the 
coupling strength was a priori unknown. They dynam- 
ically tracked and estimated the time-varying coupling 
strength for compensating against environmental fluctu- 
ation to keep it within the critical limit of stable synchro- 
nization. The critical coupling was estimated using the 
master stability function [3] which was based on assump- 
tion of identical systems. But a parameter mismatch is 
unaviodable in practical systems. A method of intermit- 
tent linear feedback coupling is found [19] to take care of 
the instability due to mismatch in delay systems, but the 
coupling did not consider the presence of a conduction de- 
lay. A delay dynamical system with conduction delay is 
of practical utility, particularly, in digital communication 
via optical-link [12]. Therefore, a unified approach for a 
robust design of delay coupling in the sense of stability 
under parameter mismatch and other external perturba- 
tion is necessary to enhance the practical use of delay 
systems. Here, we tried to address this issue starting 
with a known delay system but the coupling is a priori 
unknown. We propose a general framework of designing 
unidirectional delay coupling that can maintain stability 
of synchronization under parameter perturbation and, 
the stability is independent of the system parameters. 
We target several synchronization scenarios, complete 
synchronization (CS), antisynchronization (AS), lag syn- 
chronization (LS), anti-lag synchronization (anti-LS) and 
generalized synchronization (GS) by simply controlling 
an external parameter conveniently inserted in the defi- 
nition of the coupling. An open-plus-closed-loop (OPCL) 
coupling scheme [20] was used recently that can engi- 
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neer different desired synchronization states, CS, AS, LS 
and GS in nondelay systems. This coupling scheme has 
an advantage that it overrules many of the restrictions 
valid for other conventional coupling configurations: an- 
tiphase is observed in inversion symmetric systems [21] 
only, LS in mismatched systems [22], GS in nonidenti- 
cal or largely mismatched systems [23]. We utilize these 
common benefits of the OPCL scheme and extend it to 
delay systems in presence of coupling delay. We derive 
the stability condition of synchronization based on the 
Krasovskii-Lyapunov function theory supported by the 
Hurwitz matrix criterion [20] when the stability of syn- 
chronization becomes independent of the system parame- 
ters and is valid for mismatched oscillators too. Stability 
being free from the system parameters makes the control 
of synchronization more viable. By control of synchro- 
nization, we mean a strategy to change one to the other 
form of synchronization, either smoothly or abruptly, 
without loss of synchronization during a transition. This 
control of synchronization is really implemented by in- 
serting a scaling parameter in the design of the coupling. 
Using this parameter, a scaling of the size of a driver at- 
tractor (amplification or attenuation) is also possible at 
a response system. Additionally, we can target a type 
of mixed synchronization (MS) [20] in a response oscilla- 
tor when different pairs of state variables of the systems 
attain separate form of synchronization, CS and AS. A 
partial amplitude death [24] can also be induced in one of 
the state variables of the response system. Interesting to 
note that an amplitude death (AD) [25] can be induced 
in the response system by stabilizing its equilibrium ori- 
gin. We support the theory of the delay coupling design 
with numerical simulations of Mackey-Glass system [26] 
and a time-delayed Rossler system [27]. 

The rest of the paper is organized as follows: the gen- 
eral theory of the OPCL delay coupling in time-delayed 
systems is discussed in Section II. The stability condi- 
tion of the desired state of synchrony is derived for mis- 
matched delay systems using Krasovskii-Lyapunov func- 
tion theory and Hurwitz matrix criterion. In Section 
III, numerical examples of the Mackey-Glass system and 
the time-delayed Rossler system are described for both 
instantaneous and delay coupling. In Section IV, the 
targeting of the generalized synchronization state is dis- 
cussed with the example of time-delayed Rossler system. 
Results are concluded in Section V. 



{o-ij)nxn is a constant matrix called as the scaling matrix. 
The response system after coupling is, 

x = f(x,x T ) + D(x,x T ,y, y T ) (2) 

where D is a coupling term to be defined or designed. 
The error signal of the coupled system is defined by e = 
x — ay Tc where r c > is the coupling delay. The f(x,x T ) 
can be written as 

f(x, x T ) = f(ay Tc + e, ay T+Tr + e T ) = f(ay Tc ,ay T+Tc ) 



df 

+ — re 



df 



d{ay Tc ) d(ay T+Tc ) 



e T + ■ 



(3) 



Keeping the first order terms in (3), the error dynamics 
is approximated, 

e = x-ay Tc 



q j d f 
f(ay Tc ,ay T+Tc ) + ^— ^ e+-^— -e T + D - ay T 



d{ay Tc ) d(ay T+Tc ) 



= H 1 e+H 2 e T 
where 
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D = ay Tc - f(ay Tc ,ay T+Tc ) + {Hi 



+(H 2 - S7 — J e T 



df 



d{ay Tc ) 



d(ay T+Tc )' 



(5) 



is now defined as the OPCL delay coupling and H\ 1 H2 
are two arbitrary constant matrices. 

To derive the stability condition of synchronization, we 
first apply the Krasovskii-Lyapunov function [17], 

1 f° 

V(t) = -e T e + j e T (9)e(9)d9 (6) 

2 J - T 

where 7 is a constant. Taking time derivative of (6) and 
using the error equation (4), 

V = i[( e T i/ 1 T + e ^ 2 T )e+(ff ie +i/ 2eT ) e T ]+7(e T e-e^e r ) 



II. GENERAL THEORY OF OPCL DELAY 
COUPLING 

We consider a drive-response system where the cou- 
pling is undircctional. The driver is, 



y = f(y,y T ) + Af(y,y T ), y,y T e R n 



(1) 



where Af(y,y T ) is the mismatch parameter and y T = 
y(t—r). It drives a response system x = f(x, x T ), x,x T s 
R n , to induce a desired goal dynamics g(t) = ay(t); a = 



= E T QE 
where 
E ( " 



\{Hl + H l )+ 1 I n ^ n \H 2 



, Q 



\{H? + H 1 )+ 1 I nxn \H 2 
- 7 J n> 



(7) 
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Now we provide a proof that V < if Q < 0: 

Let Q and R be two symmetric matrices and a matrix 
S has a suitable dimension, then 



Q 
s T 



< 



if and only if both R < and Q - SR- 1 ^ < 0. 
From the above results we see that V < if 



H = 2(H 1 +H^)+4rfI n 



1 H 2 H 2 T < (8) 



The asymptotic stability of synchronization is ensured 
provided the elements of the Hi and H 2 matrices are so 
appropriately chosen that the condition (8) is satisfied. 
We construct H\ and H 2 from the Jacobian matrices of 
the system, „, ^ — , and 



Of 



™»i . If the elements of 

\ay Tc ) o(ay T+Tc ) 

the Jacobian matrices contain state variables, we replace 
them by constants hi keeping other elements (constants 
or zeros) unchanged. On the other hand if there is no 
state variable in the Jacobian, we replace some of the 
zero elements by constants pi , keeping other constants as 
it is. For a 3 x 3 Jacobian, the characteristic equation of 
a derived H matrix is A 3 + ai A 2 + a 2 X + a^ = where its 
coefficients <ij(i = 1,2,3) are defined by the elements of 
the matrix. The corresponding Routh-Hurwitz (RH) cri- 
terion [20] is given by a\ > 0, > 0, a\a 2 > when the 
H matrix has eigenvalues all with negative real parts. 
In fact, the parameter values hi and pi in Hi and H 2 
are so selected that the RH criterion is satisfied and the 
condition (8) is thereby fulfilled. It ensures asymptotic 
stability of synchronization. To complete the design of 
the OPCL coupling D, the a matrix and the coupling 
delay r c remains to be decided. The stability condition 
(8) clearly depends upon the choices of h{ and pi but 
independent of the system parameters and the choice of 
a matrix and r c . This incorporates a great flexibility in 
this coupling design in realization of several synchroniza- 
tion scenarios, CS, AS, LS, anti-LS, GS and MS in mis- 
matched systems. The coupling scheme is not difficult to 
implement in electronic circuits and laser experiments. 

As in many practical delay dynamical systems, x,y € 
R , when a is a scalar parameter, one can smoothly con- 
trol synchronization using this parameter and induce flip- 
flopping between CS (a = 1) and AS (a — —1) states, 
which is of practical use in digital encoding. A scaling 
of the size (amplification or attenuation) of the driver 
attractor at the response oscillator is possible by taking 
| a |> 1 or | a |< 1. In higher dimensional systems too, 
x,y £ R n , we can realize similar kind of scenarios and 
in addition, we can observe some new synchronization 
regimes. We elaborate using numerical examples in the 
next section. 



III. TARGETING SYNCHRONIZATION: 
NUMERICAL EXAMPLES 

We first elaborate the design of delay coupling using 
the mismatched Mackey-Glass system, 
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FIG. 1: (Color online) Synchronization region in the Hi — Hi 
plane in black in (a) for Mackey-Glass system. For r c = 0, a 
pair of time series of driver (online blue) and response (online 
red) and, the synchronization error in black lines are shown, 
(b) CS for q = 1, (c) AS for a = — 1, (d) CS with amplification 
for a = 2, and (e) AS with amplification for a — —2, (f) 
response at amplitude death for a = 0. 



V = -(a + Aa)y + 



{b + Ab)y T 
1 + V™ 



x = —ax 



1 



bx T 
ITo 



D 



(9a) 



(96) 



which are chaotic for a set of parameters a = 1,6 = 
2, r = 2 and mismatches Aa = 0.2, A6 = 0.3. The re- 
sponse system (9b) after adding the OPCL delay cou- 
pling (5) is 



-aAay Tc 



a(b + Ab)y 



T+T c 



bay T 



1 



V 10 

i/T + 7 



1 + a 1 V% 



1 - 9a 1 V+r 

+(H 1 +a)(x-ay Tc )+[H 2 -b 1 ^ T ° ](x T -ay T+Tc ) 

y i -|- ot. y r j^ Tc ) 

(10) 

For this example, Hi and H 2 are scalars, the condition 
for asymptotic stability is derived from (8), 



H 

Hi + 1 + -4- < 



(11) 



The stability region for synchronization in the H± — H 2 
plane is numerically determined as shown in the black 
region in Fig. 1(a). We then make a choice of Hi = 
—2,H 2 = —1 for our simulations. We first discuss the 
case with zero coupling delay, r c = 0. Figure 1(b) shows 
that the pair of time series of the driver y(t) (online blue 
color) and the response x{t) (online red color) are iden- 
tical in amplitude and phase for a = 1 and in a CS state 
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(a) (b) 




FIG. 2: (Color online) Lag synchronization in Mackey- Glass 
system: time series of driver (blue) and response (red) for (a) 
a = 1.0, r c = 1.0, (b) a — 1.0, r c = 3.0. Similarity measure 5* 
has a global minimum at (c) t s q = 1.0, (d) t s q = 3.0. 



since the synchronization error e(t)=x(t) — y{t) is at zero 
in black line. For a = — 1, the AS state between x(t) and 
y(t) is clear in Fig. 1(c) where the synchronization error 
e(t)=a:(t) + y(t) is found at zero in black line. Figures 
1(d) and 1(e) show the CS {a > 1) and AS (\a\ > 1) 
respectively with amplification of the driver at the re- 
sponse in both the cases. Their corresponding errors at 
zero are shown in black lines. Similarly, attenuation can 
be observed by simply taking <| a |< 1 but we do not 
elaborate here. It is also interesting that an amplitude 
death (AD) can be induced in the response system by 
taking a = as shown in Fig. 1(f). The driver y(t) is os- 
cillatory while the response x(t) ceases to oscillate and is 
stable at zero. The coupling is able to stabilize the equi- 
librium origin and thereby induces AD in the response 
system. In fact, the coupling has an inherent property to 
induce AD in a response system by stabilizing the equi- 
librium origin and even create a new equilibrium at origin 
if not present in the uncoupled system and then stabilize 
it. We like to elaborate this case in the future. 

For nonzero coupling delays r c , systems (9) and (10) 
develop a LS state instead of a CS state. In a LS state, all 
the state variables of the oscillators maintain amplitude 
correlation but shifted by the desired coupling delay r e . 
Figure 2 plots the pair of time series of the response x(t) 
(online red color) and the driver y(t) (online blue color) 
in LS for two different choices of time lag t c = 1.0,3.0. 
Figures 2(a) and 2(b) in upper panels show LS with (a) 
t c = 1.0, and (b) r c = 3.0. To confirm LS, we calculate 
a similarity measure (S)[l, 19, 21] using the pair of time 
series of y{t) and x(t), 

s _ < W) - y(t - Ts)} 2 > 

V< x?(t) >< yHt) > 



Figure 2(c) shows the similarity measure S for any delay 
t s . S has a global minimum at zero for a t s —t s q = 1.0 
which estimates the lag time between the pair of x(t) and 
y(t) in Fig. 2(a) as r c = 1.0 = r s o(dimensionless). A sim- 
ilar result is shown in Fig. 2(d) to confirm the LS with a 
larger time lag r c = 3.0 = t s q in Fig. 2(b). Thus it is pos- 
sible to induce any arbitrarily desired lag time t c in the 
drive-response system. We choose r c values one larger 
and another lower than the intrinsic time lag r and, LS 
is not affected by this choice. The efficacy of a coupling 
design for larger coupling delay is, particularly, impor- 
tant in regards to synchronization in distant regions of 
the brain with large conduction delays [28]. This knowl- 
edge can be utilized in practical situations. 

As a second example, we consider a time-delayed 
Rossler system [27]. The driving system is 

h = -Vi - e a3(t " r) (12a) 
2/2 = 2/1 + ay 2 (126) 

y-i ^ yi + {b + M)e~ v ^- T) - c (12c) 

where A6 is a mismatch. The driver system is chaotic 
for a = 0.4,6 = 2.0, c = 4.0 and < r < 0.35. In our 
simulation we take r = 0.1 and Ab = 0.05. 

After adding the OPCL coupling (5), the response sys- 
tem becomes, 

x x = -x 2 - e^- T > - a u {y 2 {t - r c ) + e ^(*-r-r„) } 
+a 22 y 2 (t - Te ) + e^a (t-T-T C ) +pi{x3 _ a33y3 ( t _ Tc) } 
+{hi+e a33y3( - t - T ^ ) }{x 3 (t~T)-a 33 y 3 (t-T~T c )} (13a) 
x 2 = a;i+ax 2 +(Q:22-aii)2/i(^-Tc) 

(136) 

x 3 = £i+6e~ X3( *" T) -a 3 3C+(a33-aii)yi(t-T c ) 

+a 33 (6 + Ab)e- y3{t ~ T - T ^ - be- 013 ™ 3 ^^ 

+{h2 + be- a33 y^- T -^}{x 3 (t - t) - a 33 y 3 (t - r - r c )} 

(13c) 

where we choose the constant matrix a — 
[an ; a 22 ; a 33 ] as a simple 
diagonal matrix. In this case, x,y G i? 3 , and hence 
Hi and H 2 matrices are 3x3 matrices. The Jacobian 
matrices of the system (12) are 

qf /0 -1 0\ ~t /0 - e M(t-r) \ 

?f= 1 a ; = 00 

dy \ 1 / ^ \ -be-y^-) J 
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FIG. 3: Synchronization region in hi — pi plane in dark black 
for delay Rossler system, hi — hi 




FIG. 4: (Color online) Mixed synchronization in delay Rossler 
system: CS between yi and Xi for an = 1 in (a) confirmed by 
yi vs. xi plot in (b), AS between yi and X2 for Q22 = — I in 
(c) confirmed by 1/2 vs. X2 plot in (d) and, (e) chaotic driver 
2/3 and corresponding response £3 at amplitude death in (f) 
for 033 = 0. 



We construct the matrices Hi and Hi from the above 
Jacobian matrices by following the rules stated above: 
Hi = [0 - 1 pi ; 1 a ; 1 0] and 
H 2 = [0 fti ;0 ; h 2 }. By a suitable 
choice of the parameters pi — —1 and hi = —2, h 2 = —2, 
the eigenvalues of the H matrix are determined as -12, 
-4, -2.4 which satisfies the Hurwitz condition(S) and en- 
sures asymptotic stability of synchronization. This choice 
is certainly not unique, other choices of parameters are 
available for which H satisfies (8) and can be easily de- 
termined as shown in the first example. The parame- 



ters in the hi — pi plane for stable synchronization are 
obtained, using numerical simulation, as shown in the 
dark region in Fig. 3 where we consider h 2 = hi. By 
this equality, we find a symmetric non-synchronization 
region in white. However, other choices of parameters 
are always possible for hi ^ h%. Clearly, a wide choices 
of pi — hi parameters is available. Once the asymp- 
totic stability of synchronization is thereby established, 
we can easily target several synchronization scenarios, 
similar to the first example, by making different choices 
of the elements of the matrix a and coupling delay r c . If 
otu = a 22 = = 1.0, with no coupling delay (r c = 0), 
systems (12) and (13) are in CS. It means all similar pairs 
of state variables of the driver and the response maintain 
a CS state. Similarly, for an = a 22 — 033 = —1.0, all 
the pairs of state variables of the systems are in AS state. 
For a numerical example, we make a different choice of 
an = 1,0^22 — ~l7 a 33 = without a coupling delay 
to reveal an uncommon synchronization scenario. As a 
result, we find that different pairs of state variables of 
the driver and the response develop different types of co- 
herence: yi(t) and Xi(t) are in CS in Fig. 4(a), y 2 (t) and 
x 2 {t) in AS in Fig.4(c). The y(t) vs. x(t) plots in Fig.4(b) 
and 4(d) confirms CS, AS respectively in similar pairs of 
state variables. The third variable 2/3 (t) of the driver is in 
chaotic state in Fig. 4(e) while the corresponding response 
variable £3 (t) is in a steady state or at partial amplitude 
death [24] in Fig. 4(f). This typical scenario is defined 
as a MS state where three different pairs of state vari- 
ables of the systems (12) and (13) exist simultaneously 
in different coherent states. Note that the system (12) 
is not inversion symmetric, even then AS is realized in 
this system. The coupling thus overrules the restriction 
of inversion symmetry as valid for diffusive coupling. 

In presence of a coupling delay, systems (12) and (13) 
will obviously develop a LS state. Instead, as a numer- 
ical example, we present a more general scenario of MS 
state with LS for a coupling delay r c = 1.0. Let us 
choose an = 1.0, a 22 = —1.0 and CV33 = 1.0 : three 
different pairs of state variables simultaneously develop 
LS, anti-LS and LS but with a common delay r c = 1.0. 
In Fig. 5(a), the time series of yi(t) (online blue color) 
and Xi{i) (online red color) are shown in LS. The corre- 
sponding xi(t) vs. yi(t — T c ) plot confirms LS in Fig. 5(b). 
Figure 5(c) shows anti-LS in the pair of time series for 
a 22 = —1.0 and a coupling delay r c = 1.0, which is con- 
firmed by the x 2 {t) vs. y 2 (t — r c ) plot in Fig. 5(d). In 
case of anti-LS, a pair of state variables are in opposite 
phase and amplitude correlated but shifted by the desired 
time lag r c . The Xs(t) and ys(t) in Fig. 5(e) also shows 
LS since we choose 0133 = 1.0. Figure 5(f) plots Xa(t) 
vs. 2/3 (t — t c ) to confirm LS between 2/3 (i) and xs(t). 
A scaling of the size of the driver attractor at response 
is also possible by making the elements of the a matrix 
larger/lower than unity. It is demonstrated that a choice 
of the elements of the a matrix is so flexible that one 
can easily target different desired synchronization sce- 
narios. The choice of the a matrix does not interfere 
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FIG. 5: (Color online) Lag synchronization in delay Rossler 
system: (a) LS for an = 1, corresponding xi(t) vs. yi(t — 
t c ) plot in (b), (c) anti-lag synchronization for 022 = — 1, 
corresponding X2(t) vs. 2/2 (£ — t c ) plot in (dO, (e) LS for 
a 3:i = 1, and xa(t) vs. yz(t — r c ) plot in (f). 




FIG. 6: (Color online) Generalized synchronization in delay 
Rossler system: (a) time series of driver y\ (blue color) and 
response xi (red color), (b) shows y\ vs X\ and (c) trans- 
form driver g\ and response variable x\ in one-to-one relation. 
Similar plots of the pair (7/2, X2), (2/2 VS.X2), and (72 vs X2 in 
(d)-(f). For the third response gs, all similar plots are shown 

m (g)-(0 



in the stability of synchronization. By varying one or 
more elements of this a matrix, a practical scheme can 
be implemented for smooth control of synchronization 
from one regime to another without loss of synchroniza- 
tion during this transition. One can even externally force 
an abrupt change in one or more elements of the a matrix 
to induce a flip-flopping of synchronization which can be 
utilized for phase-shift keying [29] in digital communica- 
tion. Presently, we use this property of the a matrix to 
demonstrate a GS scenario in the next section. 



IV. TARGETING GENERALIZED 
SYNCHRONIZATION 

A GS scenario was reported in nonidentical systems or 
largely mismatched systems under linear diffusive cou- 
pling [23]. In a GS state the response system develop a 
functional relationship with the driver. GS in delay sys- 
tems under delay coupling is also investigated [24], how- 
ever, the coupling is assumed a priori known. Instead, 
we follow the design strategy of the OPCL delay coupling 
for engineering a GS state even in slightly mismatched or 
identical delay systems. We appropriately choose the el- 
ements of the a matrix for targeting a desired functional 
relationship between the driver and the response. As a 



simple example, we choose some of the off-diagonal ele- 
ments of the scaling matrix a as constants in addition 
to the principal diagonal elements, the driver-response 
system enter into an interesting GS regime. For numer- 
ical demonstration, we use the coupled delayed Rossler 
systems (12) and (13) with an arbitrary choice of the 
scaling matrix: a = [1 — 3 — 2 ; 1 2 ; 0.5 1]. 
The goal dynamics g(t) = [gi(t) g 2 (t) 53(i)] T is de- 
fined by [ gi (t) = yi {t) - 3y 2 {t) - 2y 3 {t),g 2 (t) = y 2 (t) + 
2y 3 {t) , g 3 {t) = 0.5y 2 (t) + y3 (t) - The state variables of 
the goal dynamics or the response variables are linear 
combinations of the driver variables since the elements 
of the a matrix are chosen constant or zero. Figure 6(a) 
shows time series of the driver y\ (t) (online blue color) 
and the response X\(t) (online red color) which appar- 
ently shows no correlation between them by visual check. 
The y\{t) vs. x±(t) plot appeared to show no correlation 
in Fig. 6(b). However, the x\(t) vs. <7i(t) plot in Fig. 6(c) 
clearly shows 1:1 correlation and confirms that the tar- 
geted GS state is actually realized at the response system. 
Other pairs of state variables similarly show their corre- 
sponding targeted GS states defined by g 2 (t) and g 3 {t) in 
Fig.6(d)-(f) and Fig.6(g)-(i) respectively. Since the scal- 
ing matrix is only involved in the definition of a GS state 
and it does not disturb the stability of synchronization, 
in general, one can easily make a wide variety of choices 



7 



of the scaling matrix to define a desired functional rela- 
tionship between the driver and response. A nonlinear 
functional relationship can also be adopted to target a 
more complex GS scenario by inserting state variables of 
the driver or a third dynamical system into the a matrix 
as shown, in detail, in case of instantaneous OPCL cou- 
pling for GS [30] . The proposed design scheme is thereby 
able to engineer a varieties of desired GS state in delay 
system too in presence of delay coupling. 

V. CONCLUSION 

In conclusion, we started with a known delay system 
and a priori unknown coupling and, then proposed a 
design of unidirectional delay coupling for engineering 
various synchronization scenarios. We made necessary 
modifications in the theory of the instantaneous OPCL 
coupling for implementing in delay systems with a con- 
duction delay. We improved the stability condition of 
synchronization by supporting the Krasovskii-Lyapunov 
function theory with a Hurwitz matrix criterion and 
thereby we achieve an important benefit over the con- 
ventional methods that the stability of synchronization 
became independent of the system parameters. Another 
immediate benefit is that the OPCL delay coupling al- 
lows a flexible control of the synchronization. In fact, for 
the purpose of controlling synchronization, a scaling ma- 
trix is introduced in the coupling definition. The choice 
of the elements of the scaling matrix does not interfere 
in the stability of synchronization. One can thus switch 
or smoothly change from one regime to another with- 
out loss of stability during a transition. This is a spe- 



cial property [20, 30] of the OPCL coupling which can 
be used for practical applications. The scaling matrix is 
also used for enlarging or attenuating a driver attractor 
at a response system. This method has added advan- 
tages over the conventional methods of synchronization 
in delay system under delay coupling: it can engineer 
a desired synchronization regime in the chaotic oscilla- 
tors without restrictions of the nature of the system and 
in presence of parameter mismatch. An AS state can 
also be induced in systems without inversion symme- 
try as usually found restricted for other coupling con- 
figurations. Further, contrary to conventional coupling 
schemes, a variety of GS regimes (functional relationship) 
can be engineered in slightly mismatched oscillators by a 
wider range of choices in the elements of the scaling ma- 
trix. GS can be induced even in identical oscillators. We 
supported the theory with numerical simulations of the 
Mackey-Glass system and a delayed Rossler model. The 
OPCL scheme is thus successfully extended to delay sys- 
tems with delay coupling for targeting any of the desired 
synchronization states, CS, AS, LS, anti-LS, AD and a 
MS state. This delayed version of the OPCL coupling is 
physically implementable in time-delayed optoelectronic 
feedback loops and electronic circuits which is our future 
objective. A delay dynamical system represents many 
natural systems and since it has practical utility in laser 
based communication, a robust design of coupling, in the 
sense of stability to internal and external perturbations, 
such as the OPCL delay coupling is always important in 
search of future applications. 

S.K.D. acknowledges partial support by the BRNS, In- 
dia under grant 2009/34/26/BRNS. 
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